Acoustic surface plasmons in the noble metals Cu, Ag, and Au 
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We have performed self-consistent calculations of the dynamical response of the (111) surface of 
the noble metals Cu, Ag, and Au. Our results indicate that the partially occupied surface-state 
band in these materials yields the existence of acoustic surface plasmons with linear dispersion at 
small wave vectors. Here we demonstrate that the sound velocity of these low-energy collective 
excitations, which had already been predicted to exist in the case of Be(0001), is dictated not only 
by the Fermi velocity of the two-dimensional surface-state band but also by the nature of the decay 
and penetration of the surface-state orbitals into the solid. Our linewidth calculations indicate that 
acoustic surface plasmons should be well defined in the energy range from zero to ~ 400 meV. 

PACS numbers: 71.45.Gm, 73.20.At 



I. INTRODUCTION 

During the last decades a variety of metal surfaces, 
such as Be(0001) and the (111) surfaces of the noble met- 
als Cu, Ag, and Au, have become a testing ground for 
many experimental and theoretical investigationsi2*2i4 
These surfaces are known to support a partially occu- 
pied band of Shockley surface states with energies near 
the Fermi level. Since these states are strongly localized 
near the surface and disperse with momentum parallel to 
the surface, they can be considered to form a quasi two- 
dimensional (2D) surface-state band with a 2D Fermi en- 
ergy e|P equal to the surface-state binding energy at the 
r point. 

In the absence of the three-dimensional (3D) substrate, 
partially occupied Shockley surface states would support 
a 2D collective oscillation, the energy of this plasmon 
being given by (unless stated otherwise, atomic units are 
used, i.e., e 2 = Ti = m e = 

^2D = ^27r n 2 Dq/m 2 D, (1) 

where ri2D represents the density of occupied surface 
states: n 2 D = /tt, q represents the magnitude of 
a 2D wave vector, and m2D is a 2D effective mass. 
Eq. 0) shows that at very long wavelengths plasmons 
in a 2D electron gas have low energies; however, they 
do not affect electron-hole (e-h) and phonon dynamics 
near the Fermi level, due to their square-root dependence 
on the wave vector. Much more effective than ordinary 
2D plasmons in mediating, e.g., superconductivity would 
be the so-called acoustic plasmons with sound-like long- 
wavelength dispersion^ 

Recently, it has been demonstrated that in the pres- 
ence of the 3D substrate the dynamical screening at the 
surface provides a mechanism for the existence of a new 
acoustic collective mode, whose energy exhibits a linear 
dependence on the 2D wave vector*^ We refer to this 
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FIG. 1: Simplified model in which surface-state electrons com- 
prise a 2D sheet of interacting free electrons at z = z&. All 
other states of the semi-infinite metal are assumed to com- 
prise a plane-bounded 3D electron gas at z < 0. The metal 
surface is located at z = 0. 

mode as acoustic surface plasmon (ASP), to distinguish 
it from the conventional surface plasmon predicted by 
RitchiciiS The energy of this latter plasmon is known to 
be lu s = lo p /\/2, where uj p is the plasmon energy of a 
homogeneous electron gas of density no: lu p = (Aimo) 1 ' 2 . 

In a simplified model in which surface-state electrons 
comprise a 2D electron gas at z = z d (see Fig. QJ, while 
all other states of the semi-infinite metal comprise a 3D 
substrate at z < 0, one finds that both e-h and collec- 
tive excitations occurring within the 2D gas can be de- 
scribed with the use of an effective 2D dielectric function, 
which in the random-phase approximation (RPA) takes 
the forrn^ 

e eff(q,ui) = 1 - W(z d ,z d ;q,u;)x2D(q^), (2) 

W(z,z';q, u) being the 2D Fourier transform of the so- 
called screened interaction in the presence of the 3D sub- 
strate alonejii and XzoCff'^) being the noninteracting 
density-response function of a 2D electron gas£ 
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In the absence of the 3D substrate, W(z, z'\ q, u>) yields 
the 2D Fourier transform of the bare Coulomb interac- 
tion and e e ff(q,oj) coincides, therefore, with the RPA 
dielectric function of a 2D electron gas, which in the 
long-wavelength (q — > 0) limit has one single zero cor- 
responding to collective oscillations at u> = LJ2D ■ 

In the presence of a 3D substrate, the long-wavelength 
limit of e e ff(q,u}) has two zeros^ One zero corresponds 
to a high-frequency oscillation of energy uj 2 = uj 2 + u)\ D 
in which 2D and 3D electrons oscillate in phase with one 
another. The other zero corresponds to a low- frequency 
acoustic oscillation in which 2D and 3D electrons oscillate 
out of phase. The energy of this low-frequency mode is 
found to be of the form^ 



a vjP g, 



where vjP represents the 2D Fermi velocity 



v% D = ^2e 2 F D /m 2D , 



and 



with 



7T[7T + 2 /(%)]' 



lim W(z dl z d ; q, av 2 F D q). 

q->0 



(3) 



(4) 



(5) 



(6) 



The coefficient a (whose value depends on the electron 
density of the 3D substrate and increases with z d ) ranges 
from a constant value (on the order of 1.3 — 1.6 for metal- 
lic densities) for a 2D sheet far inside the 3D susbtrate to 
the asymptotic value \/2z d (see Ref. [T^l for a 2D sheet 
far outside the metal surface. 

In this paper, we extend the self-consistent calculations 
of the dynamical response of Be(001) reported in Ref. 
to the case of the (111) surface of the noble metals Cu, 
Ag, and Au. Our results indicate that the partially oc- 
cupied surface-state band in these materials yields the 
existence of acoustic surface plasmons whose energy is 
of the form of Eq. ©, but with an a coefficient that is 
much closer to unity than expected from the simplified 
model described above. Furthermore, we demonstrate 
that the sound velocity (v s = avjP) of this low-energy 
collective excitation is dictated not only by the Fermi ve- 
locity of the 2D surface-state band but also by the nature 
of the decay and penetration of the surface-state orbitals 
into the solid. We also investigate the width of the cor- 
responding plasmon peak, which dictates the lifetime of 
this collective excitation. 



II. THEORY 

In order to achieve a full description of the dynami- 
cal response of real metal surfaces, we first consider the 



one-dimensional potential of Ref. This allows us to 
assume translational invariance in the plane of the sur- 
face, which yields, within linear-respone theory, the fol- 
lowing expression for the electron density induced by an 
external perturbation (£> ext (z; q,u>): 

5n(z;q,u) = f dz' x(z, z'; q, u) 4> ext {z'; q, w), (7) 

x{z, z'\q, to) representing the 2D Fourier transform of 
the density-response function of our interacting many- 
electron system. 

The collective oscillations created by an external po- 
tential of the formic 



<f> e * t {z;q,«)) = -(2*/q)e* 



(8) 



can be traced to the peaks of the imaginary part of the 
so-called surface response function g(q, oj)>I£iiI 

g(q, w) = ~ J dz J dz' e« {z+z>) x(z, z'; q, u), (9) 

which at q = exhibits a pole at the conventional surface 
plasmon uj s i& 



A. Single-particle states 

The starting point of our calculations is a set of single- 
particle states V'k.raO") and energies E^ n of the form 



^k,n(r) 



■ r "i(z) 



and 



fc 2 
2m n 



(10) 



(11) 



where r = (r||,z), A is a normalization area, and 4> n {z) 
and e n are the eigenfunctions and eigenvalues of a one- 
dimensional Schrodingcr equation of the form 



2dz 2 



Vmp(z) 



4>n{z) = e n <j>n(z), 



(12) 



Vmp being the model potential described in Ref-lbl This 
potential reproduces the key features of the surface band 
structure, which in the case of the (111) surface of the 
noble metals are the presence of a band gap at the center 
of the 2D Brillouin zone (2DBZ) and the existence of 
Shockley and image states in it. 

Alternatively, for a description of the screened inter- 
action W(z, z'\ q, u>) entering Eqs. J2J and JJ]J (which ac- 
counts for the presence of 3D bulk states alone), the wave 
functions (f> n (z) and e„ can be taken to be the eigenfunc- 
tions and eigenvalues of a jellium Kohn-Sham Hamil- 
tonian of density- functional theory (DFT)ji^ which we 
evaluate in the local-density approximation (LDA) with 
the paramctrization of Perdew and Zunger^S 
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In order to solve cither Eq. i|12[l or the jcllium Kohn- 
Sham equation of DFT, we consider a thick slab with a 
given number of atomic layers and assume that the elec- 
tron density vanishes at a distance zq from either crystal 
edge^i The one-dimensional wave functions 4> n (z) are 
then expanded in a Fourier series of the form2£ 



1=1 



'n,0 



V2 
Vd 



4.i cos 



2irl 



sin 



2ttZ 
~d' 



where the distance d is given by the equation 
d = Nd + 2z , 



(13) 



(14) 



N and do being the number of atomic layers and the 
interlayer spacing, respectively; in the case of the (111) 
surfaces of the noble metals, we take N = 81, zq = 10<io, 
and l m ax corresponding to an energy of 150 eV. Due to 
the symmetry of the model potential entering Eq. (|12|l . 
the eigenfunctions (j) n (z) are easily found to be either 
even (c~ ; = 0) or odd (c+ ; = 0). 

Figure [5] shows the energies -Ek,n that we have ob- 
tained from Eq. JlT} for Cu(lll) by solving Eq. {T^J) 
as described above and using the experimental values of 
the effective masses m n of all bulk and surface surface 
states. Electronic structures for Ag(lll) and Au(lll) 
are similar. The corresponding parameters are reported 
in Table IJ 



B. Noninteracting density-response function 

Once we have an accurate description of the single- 
particle orbitals 4> n (z) and energies e n , we evaluate 
the 2D Fourier transform x°( z , z '\ <Zi w ) of noninteracting 
electrons moving in either the model potential Vmp(z) 
or the jellium effective Kohn-Sham potential of DFT: 

x\z iZ '-q,L0) = - <Pn(z)0n'(z)Mz>n'(z') 



E 



A 

n,n 

/k,n — /k+ 



(15) 



Here, the sum over n and n' includes both occupied and 
unoccupied states, rj is a positive infinitesimal, and /k.„ 
are Fermi factors, which at zero temperature are simply 



given by the Heavisidc step function 

/k,n = 0(£_F — -Ek,n) 



(16) 



ep being the Fermi energy of the solid. In particular, the 

noninteracting density-response function of the quasi-2D 

band of occupied Shockley states in the absence of the 

3D substrate can be obtained from Eq. (|15fl by omitting 

all bulk states in the sum over n and n'. 
1 
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FIG. 2: Surface electronic structure of Cu(lll). The thick 
solid line is the energy E^ } „ = e|P + fc 2 /2m2c of the Shockley 
surface state versus the 2D momentum k, as obtained with 
the measured value of the effective mass m2D ■ The grey area 
(with the upper limit at E = eo + fc 2 /2mo) represents the 
projected bulk band structure. The values of eff, ni2D, £o, 
and mo are reported in Table U] Solid lines correspond to the 
bulk states that we have obtained by employing a slab with 
81 atomic layers to simulate the semi-infinite solid. We have 
used the following parameters in the description of the model 
potential Vmp entering Eq. GIJ: Aio = -11.805, A\ = 5.14, 
A2 = 4.4204, and (3 = 2.8508 (for a description of the model 
potential see Ref.hj). Due to the presence of the band gap, for 
optical (q = 0) transitions to occur from an occupied 3D bulk 
state to an unoccupied 2D surface state the minimum energy 
^intor i s required, which decreases as the momentum transfer 
q increases. For q larger than g mln , transitions from occupied 
(unoccupied) 3D bulk states to unoccupied(occupied) surface 
states can occur at arbitrary values of the energy transfer ui. 



Introducing the one-dimensional wave functions of 
Eq. i|13|) into Eq. (|15fl . one finds the following Fourier rep- 
resentation of the noninteracting density-response func- 
tion (see Appendix A)£& 



rt/ / \ 0+/ \ f 27T?i \ / 2nn' A ^ ^ _. , . /27rn \ . ( 2-kv! A . . 

X (z,z;q,u) = }^ ^ x„> (q, w) cos I — z I cos I — —z 1+^ 2^ ^(9, "0 am I —z I sin I — « I. (17) 

n=0 n'=0 \ / \ / n —i n i—i \ / \ / 



I 

C. Interacting density-response function interacting many-electron system is obtained by solving 

In the framework of the RPA^ the 2D Fourier trans- 
form x(zi z'\ q,to) of the density-response function of an 
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FIG. 3: Energy-loss function Img(q,u})/ui of the (111) surfaces 
of the noble metals Cu, Ag, and Au, shown by solid, dashed, 
and dashed dotted lines, respectively, versus the excitation 
energy ui, as obtained from Eq. @ for q — 0.01 and 77 = 1 
meV. The vertical solid lines are located at the energies u — 
Vp q, which would correspond to Eq. with a = 1. 



the following integral equation 

X(z,z';q,w) = x°(z, z';q,u) + J dzi J dz 2 x°{z,z 1 ; q, u) 

x v(z 1 ,z 2 ;q)x(z2,z';q,uj), (18) 

the ingredients of this equation being the 2D Fourier 
transforms X°( z , z'; q, to) and v(z,z';q) of the noninter- 
acting density-response function and the bare Coulomb 
interaction, respectively. All quantities entering Eq. I|18|) 
can be represented in the form of Eq. I|17[l. which 
yields the following matrix equation for the coefficients 

n" ,n'" 

x v n n^„{q)xtmAq,u>), (19) 



v n " ,n'" (q) being the corresponding coefficients of the bare 
Coulomb interaction v(z, z'; q). 



III. RESULTS AND DISCUSSION 

A. ASP dispersion 

Figure [3] shows the energy-loss function hng(q,uj) of 
the (111) surfaces of the noble metals Cu, Ag, and Au, 
as obtained from Eq. 10 for q = 0.01. This figure shows 
the presence of a low-energy collective excitation, whose 
energy is of the form of Eq. (|3J) but with an a coefficient 
that is close to unity. Furthermore, we have carried out 
calculations of Img(g,w) for several low values of q and 
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FIG. 4: The solid line shows the energy of a well-defined 
acoustic surface plasmon of Cu(lll), as obtained from the 
maxima of our calculated surface-loss function Img(q, uj) . The 
dashed line represents the energy of an acoustic surface plas- 
mon whose linewidth starts to be considerable due to the 
presence of interband transitions. The dashed dotted line is 
the plasmon dispersion of a 2D electron gas in the absence 
of the 3D system. The grey area indicates the region of the 
(q,aj) plane (with the upper limit at uj^ — vjP q + q 2 /2rri2D) 
where e-h pairs can be created within the 2D Shockley band 
of Cu(lll). The area below the thick solid line corresponds to 
the region of momentum space where transitions between 3D 
and 2D states cannot occur. The quantities center an< l Q mm 
are determined from the surface band structure of Fig. [5] 



we have found that this low-energy collective excitation 
is indeed an acoustic surface plasmon with linear disper- 
sion, as shown in Fig. 4 for the case of Cu(lll). 

In Fig. 2] we show the energy of the acoustic surface 
plasmon of Cu(lll) versus q (solid line), as derived from 
the maxima of our calculated lmg(q, o>), together with the 
well-defined plasmon energies that we obtain when only 
the surface-state band is considered in the evaluation of 
the noninteracting density-response function of Eq. H15|) 
(dashed dotted line) . While the plasmon energies of elec- 
trons in the isolated surface-state band nicely reproduce 
in the long-wavelength (q — > 0) limit the conventional 
2D plasmon dispersion u)2D of Eq. ([TJ), the combination 
of this surface-state band with the underlying 3D system 
yields a new distinct mode whose energy lies just above 
the upper edge oj^p = vjPq + q 2 /2m 2 D of the 2D e-h pair 
continuum, as occurs in the case of Be£ Furthermore, 
Fig. 01 shows that in the long-wavelength (q — > 0) limit 
the energy of the acoustic surface plasmon in Cu(lll) is 
of the form of Eq. © but with an a coefficient that is 
considerably closer to unity than expected from Eq. (jSJ 
(see Table EJ). This discrepancy can be originated in (i) 
the absence in the simplified model leading to Eqs. © 
and JSJ) of transitions between 2D and 3D states, and (ii) 
the nature of the decay and penetration of the surface- 
state orbitals into the solid. 

In order to investigate the origin of the small differ- 
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TABLE I: 2D Fermi energy (s^P), effective mass (rri2D), and Fermi veiocity {vjf) of the Shockley surface-state band in the 
(111) surface of the noble metals Cu, Ag, and Au. £o represents the energy of the bottom of the gap at the F point, mo 
represents the effective mass of the upper bulk states at the bottom of the gap. The minimum energy transfer w^er an d 
momentum transfer q mm are those defined in Fig. 2. Also represented in this table are the values of the parameter a that 
we have obtained from our full self-consistent calculations of the surface-response function of Eq. @ and from Eq. JJJ with 
z d « 0. 
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FIG. 5: In a simplified model in which the wave function 
of surface-state electrons decays exponentially away from the 
z — Zd plane (see Fig. with a decay constant 7, acoustic 
surface plasmons are found to exist whose energy is of the 
form of Eq. with the coefficient a of Eq. JKJ replaced by 
the 7-dependent 01 coefficient that we have presented in this 
figure in the case of Cu(lll) for various values of Zd- The hor- 
izontal solid line corresponds to the coefficient a of Eq. @ 
with Zd « 0. The horizontal dotted line represents the coef- 
ficient a derived from our full self-consistent calculation that 
treats bulk and surface states on the same footing. As Zd is 
shifted from the interior of the solid towards the vacuum, the 
coefficient a increases, in agreement with Ref. 0- 



ences between the plasmon energies obtained here and 
those expected from Eqs. © and JSJ, we have carried 
out calculations of lmg(q,u;) along the lines of the sim- 
plified model leading to Eqs. © and © (see Fig. \$ but 
with the 2D electron gas of Fig. ^ replaced by a more 
realistic quasi-2D gas of electrons described by a wave 
function that decays exponentially away from the z = z c i 
plane with a decay constant 7. We have found that an 
acoustic surface plasmon is present whose energy is in- 
deed of the form of Eq. (|3J) but with an a coefficient that 
strongly depends on the decay constant 7, as shown in 
Fig. [SJ This figure demonstrates that while in the limit 
as 7 — > 00 (where the quasi-2D electron gas is indeed an 
ideal 2D sheet) the coefficient a approaches the value ex- 
pected from Eq. © (horizontal solid line) . as 7 decreases 
the dispersion of the acoustic surface plasmon approaches 
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FIG. 6: The solid line is the wave function 4>{z) of the oc- 
cupied Shockley surface state of Cu(lll) at the F point, as 
obtained by solving Eq. 112II . The dashed and dotted lines 
show wave functions of the form <j>(z) ~ exp(— f\z — Zd\) with 
Zd = —3 a.u. and two different values of the decay constant: 
7 = 0.5 Oq -1 (dashed line) and 7 = 0.25 a^ 1 (dotted line), ao 
is the Bohr radius: ao = 0.529 A. 



(for all negative values of Zd) the more realistic situation 
where a is close to unity (horizontal dotted line). A 
comparison between the model wave functions that we 
have used in this calculation and the actual surface-state 
wave functions that are involved in the full calculation 
of Figs. 13 and 01 is presented in Fig. Although the 
model wave functions do not reproduce the actual shape 
of the surface-state wave function, a finite penetration of 
the model surface-state wave functions into the solid al- 
lows the formation of an acoustic surface plasmon whose 
sound velocity is very close to the Fermi velocity of the 
2D surface-state band (a ~ 1), as predicted by our more 
realistic calculation. Indeed, the finite penetration of the 
surface-state wave function into the solid provides a more 
complete screening of the quasi-2D collective excitations 
by the surrounding 3D substrate, which brings the acous- 
tic surface plasmon closer to the upper edge of the 2D c-h 
pair continuum (a — > 1)^ 
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FIG. 7: Solid and open symbols represent the width at half 
maximum of acoustic surface plasmons in Cu(lll) versus the 
plasmon energy, as obtained from the imaginary part of the 
surface-response function g(q,u)) of Eq. © with (solid sym- 
bols) and without (open symbols) inclusion of transitions be- 
tween 2D and 3D states in the evaluations of the noninter- 
acting density-response function of Eq. I)15ft and for various 
values of the parameter r\: 0.1, 1, and 5 meV. The solid and 
dashed lines represent fits from the solid and open symbols, re- 
spectively. For this surface, the threshold for interband tran- 
sitions between 2D and 3D states occurs at ~ 110 meV (see 
Fig.EJ. 

B. ASP linewidth 

Finally, we have carried out lifetime calculations of the 
acoustic surface plasmon, as derived from the width of 
the imaginary part of the full surface-response function 
of Eq. 10. Fig. [7| shows the results we have obtained 
for Cu(lll) with (solid symbols and lines) and without 
(open symbols and dashed lines) inclusion of transitions 
between 2D and 3D states in the evaluation of the non- 
interacting density-response function of Eq. l|15f) and for 
three different values of the parameter rj. Plasmon decay 
can occur by exciting e-h pairs either through transitions 
between 2D and 3D states, which would not be present in 
the model leading to Eqs. © and (JSJ), or through transi- 
tions within the 3D continuum of bulk states^ At small 
energies below the threshold of u> ~ 110 meV (see Fig.0J, 
where acoustic surface plasmons can only decay by excit- 
ing e-h pairs within the 3D continuum of bulk states (the 
solid and dashed lines of Fig. coincide), the linewidth 
is entirely determined by the choice of the parameter 
77, showing that at these low energies the impact of in- 
traband transitions between 3D bulk states is negligibly 
small. As the plasmon energy increases, there is a small 
contribution to the plasmon linewidth from transitions 
between 2D and 3D states (the difference between solid 
and dashed lines)^! and an increasing contribution from 
intraband 3D transitions yielding a finite linewidth which 
for 77 < 0.1 cV is not sensitive to the precise value of ?/ 
employed, but which still allows the formation of a well- 



defined acoustic-surface collective excitation for plasmon 
energies at least up to ~ 400 meV. A similar behavior is 
observed in the case of Ag and Au. 



IV. SUMMARY AND CONCLUSIONS 

We have carried out self-consistent calculations of the 
surface- loss function of the (111) surfaces of the noble 
metals Cu, Ag, and Au, by considering a one-dimensional 
potential that describes the main features of the surface 
band structure. We have found that the partially oc- 
cupied surface-state band in these materials yields the 
existence of acoustic surface plasmons, as it had already 
been demonstrated to occur in the case of Bc(0001)£ The 
energy of these collective excitations has been shown to 
exhibit linear dispersion at small wave vectors, with the 
sound velocity being very close to the Fermi velocity of 
the 2D surface-state band and considerably smaller than 
expected from a simplified model in which surface-state 
electrons comprise a 2D electron gas while all other states 
of the semi-infinite metal comprise a 3D substrate. 

The origin of the differences between the plasmon en- 
ergies obtained here and those expected from simplified 
models has been investigated by performing simplified 
calculations with the ideal 2D sheet of Shockley electrons 
replaced by a more realistic quasi-2D gas of electrons 
whose wave functions decay exponentially away from the 
surface. These calculations have been found to yield an 
acoustic surface plasmon whose energy is linear in the 
magnitude of the wave vector, the sound velocity being 
dictated not only by the Fermi velocity of the 2D surface- 
state band but also by the nature of the decay and pen- 
etration of the surface-state orbitals into the solid. With 
an appropriate choice of the exponential decay of surface- 
state wave functions these simplified calculations (which 
do not account for transitions between 2D and 3D states) 
accurately account for the energy dispersion of acoustic 
surface plasmons, which indicates that the impact of in- 
terband transitions between 2D and 3D states on the 
ASP's energy dispersion is negligible. 

We have also carried out self-consistent calculations of 
the linewidth of acoustic surface plasmons in the (111) 
surfaces of the noble metals. We have found that while 
the impact of interband transitions between 2D and 3D 
states is small, intraband transitions between 3D bulk 
states contribute considerably to the finite linewidth of 
acoustic surface plasmons, which are found to represent 
a well-defined acoustic collective excitation for plasmon 
energies at least up to ~ 400 meV. 

Finally, we note that as in the case of conventional 
surface plasmons, acoustic surface plasmons should also 
be expected to be excited by light, as discussed recently^ 
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density-response function of Eq. (|17(l , as obtained by in- 
troducing the one-dimensional wave functions of Eq. 
into Eq. (|15(1 . Replacing the sum over k in Eq. (|1 by 
an integral, we find 
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Here we give explicit expressions for the coefficients 



Xnn'(9j w ) entering the expansion of the noninteracting 
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q 2 k 2 ( q 2 k 2 k 2 

+ £/' — e; + ±lo ±irj 



2 

mi 



\2to/' 



2m;/ 2m; 



-, -1/2 



(A4) 



(A5) 



,+ „+ 



^ EE 1 



n+n'' ,n' 



J n-\-n' ,n 



n'-\-n" ,n J 



(A6) 



g: 



2 E E l 



-f-n" ,n' 



+ n +7 



+n" ,n)i 



(A7) 



and 
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